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Abstract 

The growth of a diffusion limited aggregation (DLA) cluster with mass M 

and radius of gyration R is described by a set of growth probabilities {pi}, 

where pi is the probability that the perimeter site i will be the next to grow. 

We introduce the joint distribution N(a, x, M), where N(a, x, M)dadx is the 

number of perimeter sites with a- values in the range a < ati < a + da ( "a- 

sites") and located in the annulus [x, x+dx] around the cluster seed. Here, 

a,i = — lnp,/lnii! if pi > 0, x = Ti/R, and r*j is the distance of site i from 

the seed of the DLA cluster. We use N(a, x, M) to relate multifractal and 

multiscaling properties of DLA. In particular, we find that for large M the 

location of the a-sites is peaked around a fixed value x(a); in contrast, the 

perimeter sites with Pi = are uniformly distributed over the DLA cluster. 
61.43.Hv, 68.70+w, 05.40+j, 81.10.-h 
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I. INTRODUCTION 



The growth of a diffusion limited aggregation (DLA) | l [14|] cluster with mass M is 
described by the set of growth probabilities {p^ |[T5HT^] where pi is the probability that 
perimeter site i will be the next to grow. One way to analyze the set {p{\ is by calculating 
the "growth-probability distribution" n(a, M), where n(a, M)da is the number of perimeter 
sites with a < ai < a + da, 

ai = - In pi / ha. R, (1) 



and R is the radius of gyration of the cluster |18| , |19| . We call a-sites those sites which are 



characterized by the same value of a. The main motivation for studying the distribution 
n(a,M) is its relation to the multifractal spectrum f(a) |20| , ^I| of the "measure" {pi}- We 
define f(a, M) through 

n(a, M) = M ufia ' M \ (2) 

where v is the inverse fractal dimension of DLA, R ~ M v . If for large M f(ot, M) converges 
to an M-independent function f(a), then f(a) is usually called the multifractal spectrum. 
For 2D DLA, there exist several studies p2|-f2~if proposing different convergence behaviors 
and functional forms of f{ot) in the limit M — > oo. However, these considerations will not 
be essential for our arguments concerning the relation of multiscaling and multifractality. 
Henceforth we will only assume that some f{ot) exists. 

During the process of calculating n(a,M), the information about the location of the 
a-sites is lost. However, some information about the location of the growth sites and their 
associated values pi may be obtained from the Plischke-Racz probability P(x, M) |[25|| , where 
P(x, M)dx is the probability that the next particle will be deposited at a rescaled distance 
x < Xi < x + dx. Here x\ = Ti/R and is the distance from the cluster seed. For DLA, the 



function P(x, M) displays a peak at a constant value x of the deposition radius [25 
A simple form for P(x, M) is the Gaussian, 
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where £ denotes the mean square width of the deposition zone. Plischke and Racz [25 



suggest that £ ~ M v> u ', z/ < z/, where v is the inverse fractal dimension of DLA, R ~ M v , 



and i/' an independent exponent. However, Meakin and Sander [2E[] find that v' approaches 
v a.s M increases. They also argue that v = v 1 in the limit M —>■ oo. 
If (x/£) 2 ~ 2clnM |27|] with constant c, then P(x,M) takes the form 

P(x, M) = -^M-^CprCx), (4) 

where from ([|) 0(a;) = c(x/x — l) 2 and Cpr(x) = 1. In general, if Cpr(x) is a generic 
function of x and is independent of x we have conventional scaling while if <p(x) is 
x-dependent we say that we have multiscaling. 

Next, we introduce the annular density pa(x,M), where Pa{x, M)dx is the number of 
particles in an annulus [x, x + dx] and related to the conventional particle density p(r, M) 
by pa{x, M)dx = 27rrp(r, M)dr. Since the change of pa{x,M) with increasing cluster mass 
is given by P(x, M) p5[[ , i.e., 





dM 



p A (x,M) = P(x,M), (5) 



it was suggested |28j to write pa(x, M) in the same multiscaling form as Eq. ([|), 

PA (x,M)=r D ^C p (x), (6) 

where D(x) is the fractal dimension for a thin annulus with average radius x, and C p (x) is 
an amplitude. Note that if <fi = <f>(x) in Eq. (^), then also D = D(x) in (|6]). Multiscaling in 
(H) has been supported by simulations [Q . Whereas multiscaling for clusters with M < 10 6 



has been confirmed very recently by P. Ossadnik [p^j , the same study analyzes one very 
large off- lattice cluster of M = 5 x 10 7 arriving at an ambiguous result, which is consistent 
with both multi- and standard scaling of pa(x, M). 

As demonstrated in Ref. multiscaling results if the a-sites are "localized' in space. 
Here, we study the nature of the "localization" of the a-sites and the non-localized behavior 



of the pi = ("dead") perimeter sites (Sees. ||, The consequences for the multiscaling 
hypothesis of P(x,M) and Pa(x, M) will be discussed in Sec. fVj. Moreover, we introduce 
the notion of a multifractal spectrum in an annulus and find an intriguing combination of 
multifractal and multiscaling properties (Sec. |V|) P0fl . 

II. THE JOINT DISTRIBUTION FUNCTION 

In this section we introduce the joint distribution function N(a,x,M), where 
N(a,x, M)dadx is the number of perimeter sites with p, > such that a < oij < a + da 
and within the annulus [x, x+dx]. The distribution N(a,x, M) can be related to the three 
functions discussed in the introduction: 

(a) n(a,M). By integration of N(a,x, M) over x, we have, 

n(a,M) = J dxN(a,x,M). (7) 

(b) P(x, M). If we use Eq. ([[]) together with the relation R = M v — which reflects the 
well-established fractal structure of DLA — to write the growth probability as M~ ua , then 

P(x,M) = J daN(a,x,M)M- ua . (8) 

(c) pa(x,M). One possibility to express pa(x, M) in terms of N(a,x, M) is by using 
Eqs. (d) and (j|). Integration of Eq. @ with respect to M yields 

i-M r- 

Pa(x,M) = J^ dM' J daN(a,x,M')(M')- ua (9a) 

However, we would like to point out another relationship that does not involve an integration 
over the growth history of the cluster. First, note that N(a, x, M) only describes perimeter 
sites with growth probability pi > (alive perimeter sites - - in contrast to dead sites 
with pi = 0, see Fig. |1|). We denote the annular density profile of alive perimeter sites 
by Pa\x,M) and that of dead perimeter sites by p^\x,M) } both defined in analogy to 
Pa{x,M), which describes the density of cluster sites. The sum p^\x,M) + p^\x,M) 
describes all perimeter sites of clusters of mass M. Since DLA is a treelike fractal object 
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|31| , we expect the annular density of all perimeter sites to be proportional to the annular 
density of cluster sites Pa(x, M). Furthermore, as one can see from Figs. |^ and ^, the spatial 
distributions of dead and alive sites are in good approximation proportional to each other, 
i.e., Pa(x,M) ~ Pa\x,M); a detailed discussion is given in Appendix [A[ Consequently, 
we expect 

p A (x,M) ~ f daN(a,x,M). (9b) 



III. SIMULATION APPROACH 

A calculation of the joint distribution N(a, x, M) shows that the x dependence is approx- 
imately Gaussian, centered around a value x = (x) with variance £ 2 (a,M) = (x 2 ) — {x) 2 , 
i.e., 



N(a, x, M) oc exp 

/km 



(x — x(a)Y 



(10) 



2£ 2 (a,M) 

Here, the brackets ( ) indicate an expectation value with respect to the empirical distribution 
N(a,x,M), i.e., (f(x)) = f dxf(x)N(a,x,M)/ f dxN(a,x,M). 

In Fig. § we show N(a, x, M) with M = 20 000 for a values of 1.1,1.5,1.9,2.3 and 
2.8 vs. ±(x — x) 2 , where the + sign applies if x > x and — otherwise. The Gaussian 
approximation (0) is justified, since we observe for positive and negative abscissa values 
an approximate straight line behavior of N(a,M). For each a, the modulus of the slope 
m{a) of these lines relates to the width of the Gaussian, |m(a)| = l/2£ 2 (a,M). Apparently, 
£(a, M) increases with a. In other words, highly screened growth sites are less localized 
than the exposed growth sites in the active region of the cluster that are characterized by 
small values of a. The approximation is worse for x < x, where especially for large values 
of a the presence of the cluster center distorts the pure Gaussian behavior |53|. 



In Fig. H we demonstrate that a-sites are located in approximate annuli around the 
center of the cluster by displaying all the "live" (pi > 0) perimeter sites of 18 superposed 
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off-lattice DLA clusters of M = 20 000 for three distinct values of a. In contrast, the dead 
sites of a cluster, as shown in Fig. |2|a, are distributed with a density proportional to that of 
the alive sites (cf. App. which, for comparison, are displayed in Fig. Qb. 

The M dependence of the width £(a, M) is shown in Fig. ^|a. We cannot identify an 
unequivocal limit behavior as a function of M. For certain values of a, £(a, M) seems to 
decrease as M — > oo — for others an increase is observed. The best statistics are obtained 
for a ~ 1.36, for which we find a decreasing width. Since asymptotically £(a, M) can 
certainly not increase — which would correspond to the statement that the growth zone 
would become larger than the cluster itself — we are clearly still in a mass regime where 
finite size effects play an important role. Thus, we will have to discuss several possibilities 
for the M-dependence of £(a, M) in the following sections. 

The M dependence of the mean x(a, M) is plotted in Fig. ^b. Unlike in the case of 
the width £(a,M), here the limit behavior for large M is manifest. For small values of a, 
x(a, M) decreases with M while for large a, x(a, M) increases with M. However, in both 
cases x(a, x, M) converges towards an M-independent limit x(a), as shown in Fig. []c. Note 
that x(a) is a monotonically decreasing function and thus invertible. The decrease of x(a) 
with a results from the stronger screening (large a) in the interior of the cluster (small x). 



IV. MULTISCALING OF P(x,M) AND p A (x, M) 



Next we discuss the functional form of the new distribution function N(a,x, M). As 
we have shown above, our calculations are consistent with the possibility that for large M, 



N(a, x, M) is a Gaussian |32 



m x, 



N(a,x,M) 



■■ exp 



(x — x(a)) 2 
2£ 2 («,M) 



(11) 



/27r£ 2 (a,M) 

In (P|) we write x = x(a), since we have found that x for large M depends only on a. The 
term M 1 "^^ in the prefactor of the Gaussian ensures that n(a, M) has the multifractal form 
(@), as one can see by integration with respect to x. 
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We next discuss the implications for multiscaling of several possibilities for the functional 
dependence of the width M) on a and M. Preserving the multifractal properties of the 
distribution, we explore three mutually exclusive cases for the mass dependence of £(a, M): 

(A) Constant width, £(at, M) = A(a). A constant width corresponds to the possibility 
that both the average location of the a-sites and the width of the growth zone are not mass 
dependent, implying that both length scales are proportional to the cluster radius R. 

(B) "Strong localization," £(a, M) = A(a)M~ y , y > 0, corresponds to a fast decrease of 
the width of the spatial distribution of the a-sites. Asymptotically, the a-sites are located 
at the same distance x(a) from the cluster seed. 

(C) "Weak localization," £(a, M) = A(a) / y/\n M , where the width still tends to zero, 
but in a logarithmic fashion with intriguing implications for the scaling behavior of P(x, M) 
and pa(x, M). 



A. Constant Width: f (a, M) = A(a) 



Substituting relation ([LTD into @, we obtain for the Plischke-Racz probability 

M v fWM~ va 



r 

P{x,M) = / da- 



x exp 



2nA 2 (a) 

(x — x(a)) 2 



2A 2 (a) 



(12) 



Performing a steepest descent analysis of fll2), we now calculate the value a* of a at which 
the integrand is maximal. In the M — > oo limit, the resulting condition for a* is 

d 



da 



/(«)!«. = i. 



(13) 



The value of a* which satisfies Eq. ( p^) is known to be unity ||33|| , moreover, for a = 1 we 
have f(u) = 1. Thus, in case (A), P(x,M) is M-independent and has Gaussian shape with 
constant width £ = A(l), i.e., 



P(x,M) 



1 



2vrv4 2 (i; 



exp 



(x - x(l)f 
2A 2 (1) 



(14) 



Eq. (fn|) does not explicitly depend on M and thus P(x, M) obeys standard scaling. 

We perform a similar analysis to evaluate the expression ( [9b]) for the density profile of 
the cluster, 

(x — x(a)) 2 



Pa(x,M) 



M vf{a) 
da — , exp 

2-KA 2 {a) 



2A 2 (a] 



(15) 



Here, the condition for the saddle-point value a* in the limit of large M becomes, 

= (16) 
The maximum /(a*) has the value f{a*) = 1/V, so that the density can be written as 



p A (x,M) 



M 



27cA 2 (a* 



: exp 



(x - x(a*)) 2 
2A 2 (a*) 



(17) 



Formally, Eq. ([17]) can be cast into the form of Eq. 



p A (x,M)~R 1 /» r -^-exp 



(x - x(a*)) 2 
2A 2 (a*) 



r 1/u C p (x) 



(18) 



Note that in case (A) the exponent of r is independent of x and thus also pa(x,M) obeys 
standard scaling. 

It is instructive to note the meaning of the values a* from Eqs. ( P~5| ) and flTB| ) in the 
multifractal spectrum. The probability to grow at a site with a specific value of a is maximal 
when the product of growth probability and number of sites with this probability, and thus 
f(a) — a, is maximal. As can be seen by differentiating with respect to a this condition 
is equivalent to relation (|13|), from which results the dominant value a* that controls the 
P(x,M) integral flT2|). In contrast, the mass distribution p A (x,M) (|I5D in the cluster is 
controlled by the a value corresponding to the maximum of f(a), which is the fractal 
dimension of the set comprising the "majority" of growth sites. 

However, the simple possibility that the width M) of N(a,x, M) is independent of 
M appears questionable. First, a recent calculation ]33| indicates that the width of P(x, M) 
decreases with M. Second, simulation results for off-lattice DLA clusters with M up to 
5 x 10 7 are consistent with the multiscaling relation (^) for pa{x, M) which precludes an M 
independent width £(a,M). 
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B. Strong Localization: £(a,M) = A(a)M~ y 



The implications for case (B), £(a,M) = A(a)M y , are quite different from case (A). 



Now, the analog of Eq. (12) for the Plischke-Racz probability becomes 

J^vf{a)-ua 



P(x, M) = J da 



^27iA 2 (a)M- 2 y 



x exp 



-M 2y 



[x — x(a)) 2 
2A 2 {a) 



(19) 



For sufficiently large M, the Gaussian in the integrand tends to a delta function centered at 



x, 



x = x (a . 



(20) 



Given an annulus x, there is only one value of a = a(x) given by the inverse function of 
x(a). - - the monotonicity of x{a) guarantees the existence of a(x) (cf. Sec. [TTT] and Fig. 

§• 

Using (|20"D, we can compare P(x,M), Eq. (0), to the multiscaling form (f|) and find 
that 



(x) = —uf(a(x)) + vat(x). 



(21) 



Thus, in contrast to case (A), case (B) results in multiscaling. 

In the same fashion, from ( pa]) and (0) we can demonstrate multiscaling for the density 
profile of the cluster. Again the gaussian in (|TT|) tends to a delta function centered at 
x = x(a) and thus ( p0| ) also determines a(x) for the density profile. The resulting density 
profile Pa(x, M) is 



p A {x,M) ~ M v ^ a) ~ r f{a{x)) C p {x), 



(22) 



and displays multiscaling as in Eq. (0) with 



D(x)=f(a(x)), 



(23) 



in agreement with the result of Ref [27] . 



The results ( pT|) and ( p2|) are not altogether surprising. In fact, the distribution 



N(a,x, M) for fixed a tends to a 5-function as M —>■ oo. In this limit, almost all the 
sites with a specific a- value are located at distance x(a) from the cluster seed — such that 
we refer to case (ii) as "strong localization." Vice versa, a specific location x singles out an 
a value a(x). From f{a) we then obtain the fractal dimension of the set of these a(x)-sites. 
Eq. (|22D can now be understood just as the usual relationship between mass and extension 
of a fractal object and, similarly, M~^ x \ which describes the probability of deposition at x, 
is just the product of the growth probability M~ ua ^ at x and the multiplicity of the growth 
sites at x, M uf ^\ 

In the case of strong localization it is particularly simple to obtain an understanding 
of the relationship between cluster structure and the distribution of growth probabilities. 
The large a part of f{ot) reflects the properties of the frozen region (small x) of the cluster, 
where the Pi are so small that effectively no further growth will occur. One expects that 
the mass distribution in the frozen region is characterized by the fractal dimension of DLA. 
This assumption is supported by the results of Refs. p5| , pU| , pi| . However, since a(x) is 
not constant for small x, and if our assumption ( |TTj ) for the form of the joint distribution 
function N(a,x, M) is still valid, then f(a) has to be independent of a for large a. In fact, 
the phenomenon of a "phase transition" |]35|j36|] in DLA is consistent with such a behavior 



For case (C), M) = A(a)vmM. The exponential in the integrands of (|T2| ) and ([19]) 
turns into a power-law, 



of f(a) @. 



C. Weak Localization: £(a,M) = A(a)/y/]n.M 




(24) 



Given a value x, a steepest descent analysis of this integral yields 
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d ( , . (x — x(a)) 



2 N 



uf(a) - 



(25) 



da V ' 2A 2 (a) 
as condition for the value a* maximizing the integrand. As we vary x, the changing value 
a* defines a function ap(x) which enables us to write the Plischke-Racz probability as 

P(x, M) = Cp R (x)M-^ x) (26) 

where 

(f)(x) = -vf{a P {x)) + ua P (x) + 2 A 2 (a P (x)) ' ^ ' 

and Cpr(x) is an amplitude. By comparison of relation (|26|) to Eq. (^) we see that case (C) 
like case (B) results in multiscaling, but with a much more complex multiscaling "exponent" 
4>{x). 

Although the width £(a, M) of N(a, x, M) in case (C) still approaches zero for large M, 
we note that, in contrast to case (B), a* is no longer "characteristic" for the shell x — in 
the sense that Eq. (^) no longer holds. Thus, we refer to case (C) as a case of "weak 
localization." 



To analyze pa(x, M), we first write the condition for a*. In analogy to (p5|) 



d ( , . (x — x(a)) 



2 N 



uf(a) - 



0. (28) 



da V 2A 2 (a) 
The a* values satisfying Eq. ( |2"B] ) define a function a p (x), when x is varied. Unlike the 
strong localization case (B), a p (x) differs from ap(x). We can use the function a p (x) and 
the method developed in Eq. (|18D to express Pa{x, M) as 

p A (x,M) = C p (x)r D ^, (29) 

where C p (x) is an amplitude and the multiscaling exponent 

fl <*> = 1{a ^ - ^OT " (30) 

The numerical data presented in Sec. |TTT| are not sufficient to distinguish between cases 
(A) through (C) for all values of a. However, for the particular value a = 1.36, for which 
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the statistics is good, we find the data consistent with case (C). In the following section we 
will provide further support for the multiscaling idea. 

In the preceding discussion we have seen that the existence of two different length scales, 
namely the average location x(a) and the width M) in cases (B) and (C) leads to 
multiscaling for both P(x, M) and Pa(%, M). In contrast, we do not find multiscaling in case 
(A), where only one length scale is present — both x(a) and M) are M independent. 



V. MULTIFRACTALITY 

Another interesting quantity — possibly accessible to experimental measurements — is 
the scaling behavior of "moments" of N(a, x, M). In this section, we will define a multifractal 
analysis for the Pi contained in an annulus at distance x in analogy to the multifractal 
formalism presented, e.g., in Refs. pO| , pT| . 

Usually, a multifractal analysis is performed on a set of numbers that are normalized. 
However, the sum of the growth probabilities Pi restricted to an annulus is less than one, 
since the Pi are normalized with respect to the entire cluster. Here, we will base our analysis 
on the unnormalized set of Pi at a specific distance x from the cluster seed and postpone to 
App. |B| a discussion of what happens if we use normalized probabilities instead. 

Our first step is to define a "partition function" 

Z(q,x,M) = J daN(a,x,M)M- qua J (31) 



where q is an arbitrary real number p0| , |3~7|| . The function Z(q, x } M) can also be considered 
the q th "moment" of the distribution N(a,x, M). Second, we define the scaling indices 
r(q,x) as a function of q for different x by 

T(q,x)= hm ln Z [ q t M) - (32) 

If r(q, x) is a linear function of q, then conventional "gap scaling" is obtained, while otherwise 
we call the measure underlying the "moments" Z(q, x, M) multifractal. It is then convenient 
to introduce the Legendre transform fi{a,x) of r(q,x) 
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d 

h(a L ,x) = qa L - r(q,x), where a L = —r(g, x). (33) 
The quantity /l(q;l) can be interpreted as the fractal dimension of the set of points 



characterized by ai |20| in the annulus x. We analyze the scaling behavior of the moments 
by performing a steepest descent analysis of (|3T]). Substituting N(a,x, M) in its analytical 
form ( |TT| ) into Eq. ( |3"T| ) yields 

Jy[vf(a)-vqa 



da . exp 
•/27r^ 2 (a,M) 



(x-x(q)) 2 ! 

2e 2 («,M)J- 1 J 

As in Sec. [TV] we will now discuss the consequences of different asymptotic behavior of the 
width of the growth zone of the cluster. 

A. Constant Width: f (a, M) = A(a) 

For the case of constant width (cf. Sec. |IV A| ), we find that the dominant contribution 
to the integral (|4j) arises at a value a* = a* (I) given by 

(35) 

From Eq. (|35D , we can anticipate the result, 

f L (oi L ,x) = f(a L ), (36) 

which we obtain after consideration of Eqs. (|32|) and (^). Eq. ([36]) states that /or case (A) 
the full, unaltered multifractal spectrum is found in all annuli x. Thus, given the validity of 
Eq. flTTD, a constant width of the growth zone can only be maintained, if sites with both 



low and high growth probabilities are distributed evenly in the cluster . 

However, due to the the screening of the interior regions of the cluster, the growth 
probabilities for small x are significantly smaller than at the exposed sites on the exterior of 
the cluster. These smaller probabilities result in a shift of the distributions N(a, x, M) to 
larger values of a for small x. Thus, the multifractal spectrum /l(«l) differs for different x 
(cf. [|37| for the 3D case), in contrast to the result (|36|) above. 
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B. Strong Localization: f (a, M) = A(a)M~ y 



For the strong localization case (Sec. |1V B|) , the Gaussian term in ( |34"D tends to a 5- 
function localized at x given by Eq. (|20| ) whose inversion gives a = a(x). Now, we use for 
Z(q, x, M) in Eq. (^) the value of the integrand in ([31]) at a — a(x). The resulting scaling 
indices r(q, x), 

r(q, x) = qa(x) - f(a(x)), (37) 
are linearly dependent on q for fixed x. Thus, in case (B) gap scaling of the moments 



Z(q,x,M) results. The /l(«l) corresponding to (|37|) is the point [a L = a(x); /l(«) = 
/(«(*)) ]• 

The absence of multifractality in a given annulus x in case (B) is not surprising, since 
the strong localization of a-sites implies a very narrow spread of a values within a specific 
shell x. Narrow distributions, however, typically display gap scaling of their moments. 

C. Weak Localization: £(a, M) = A(a)/VlnM 

Again, the weak localization case is quite different. The value a* maximizing the inte- 
grand of fl3~Tl) is given by 



d ( , . (x — x(a)) 



2 s 



uf{a) - 



vq. (38) 



da V 2A 2 (a) 
In contrast to Eq. (|35D, where the solution depends on q only, here a* is a function of both 
q and x, and we write a* = a q (x). From Eq. fl32|), we obtain 

r(g,x) = vf[a q (x)) - va q [x)q 2 A 2 (a (g)) ' ^ ^ 



In (|39| ) r(g, x) is a nonlinear function of both q and x, so that its Legendre transform ( |33| ) 
depends on x. Thus, in case (C), the multifractal spectrum fL(ai,x) is a function of the 
location of the annulus x. In contrast to case (B), the width £(a,M) approaches zero so 
slow that the multifractality of case (A) is not destroyed as in case (B), but altered in its 
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character. In fact, the x dependence of the multifractal spectrum fi,(a,x) is the hallmark 
of multiscaling as encountered for P(x,M) and pa(x,M) in case (C), Sec. [TV} 

In Fig. [7] we display r(q, x) for different values of x as a function of q ||35|| . 

For large \q\ the function r(q,x) tends to straight lines with different slopes. The def- 
initions ([H],[J^) show that for q < the slope is determined by the mass dependence of 
the smallest growth probability p m \ n {x) within the annulus x and for q > by the mass 
dependence of the largest growth probability p ma x(^), respectively. Especially in the region 
< q < 1, we observe a pronounced curvature of r(q,x). Since the strong localization case 
predicts a linear behavior of r(q, x) over the entire range of q values, we conclude that our 
findings disfavor strong localization. 

Moreover, we see that r(q,x) displays variation with the parameter x. If x becomes 
smaller, both p m i n (x) and p max (x) as functions of the cluster mass decay faster, because the 
interior frozen regions of the cluster are screened stronger. Thus, the two asymptotic slopes 
(for q — > ±oo) of r(q, x) increase. For example, for any given q, r(q, x = 1.9) has everywhere 
a slope less than r(q,x = 0.7). As a consequence, also the Legendre transforms of r(q,x) 
are x dependent. This finding is consistent with the weak localization case (C). However, 
due to the comparatively small clusters that we have analyzed, we cannot exclude cases (A) 
or (B). 

In Fig. [5] we see that r(q, x) only starts to change appreciably for quite large x-values 
around x = 1.5. A similar phenomenon is observed in the numerical multiscaling analysis 
of the annular density of the cluster. There, the function D(x) is approximately constant 
= \jv up to x values of similar magnitude before D(x) drops to zero over comparatively 
small range. 

For 3D off-lattice DLA, the multifractal properties of the growth probabilities {p{\ in an 
annulus were calculated in ||37|| . Although no statement about N(a,x, M) for 3D DLA was 
made, the results display qualitatively the behavior predicted above in the cases (A) and 
(C) for 2D DLA. 
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VI. CONCLUSION 



We have discussed different possibilities for the analytical form of the joint distribution 
function N(a,x, M). For cases (B) and (C), where two different length scales enter into 
N(ot, x, M), we find multiscaling behavior of the Plischke-Racz probability P(x, M) and the 
density profile Pa{x, M) of DLA clusters. Moreover, we find that the scaling behavior of the 
moments of the growth probability distribution constrained to an annulus x is different in all 
cases. For (A) and (C), we encounter multifractality, which in case (C) bears an additional 
feature typical for multiscaling: the x dependence of a scaling "exponent." Our data and 
previous work is consistent with case (C), although further numerical work to clarify the 
rich scaling properties of N(a, x, M) is clearly desirable. 

Our results for all the different cases discussed above are summarized in Table |. 
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APPENDIX A: SPATIAL DISTRIBUTION OF DEAD AND ALIVE PERIMETER 

SITES 

In this appendix, we discuss the properties of the distribution of dead perimeter sites 
|19| . These sites with growth probability exactly equal to zero result because specific con- 
figurations of cluster sites enclose perimeter sites in such a way that they can no longer be 
reached from the exterior of the DLA cluster. If dead sites are predominantly formed due 
to specific local configurations, then the fraction of perimeter sites that are dead will not 
change as the cluster grows []nj. Moreover, the spatial distribution of dead sites is then 
proportional to the distribution of alive sites, i.e., 
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p { *\x,M)~p i *\x,M)~ PA ( X7 M). (Al) 

On the other hand, if the dead sites were found mainly near the region of small pi, then 
p^ (x, M) would be shifted towards the center of the cluster. In Fig. |, we notice that the 
shape of both distributions looks almost identical. The similarity in form is an evidence 
favorable to the above stated "local configuration argument," and shows that the dead sites 
are uniformly distributed over the cluster. 

The simplest local configuration producing a dead site is the "L" configuration shown in 
Fig. Qa. In general, also more complicated configurations will contribute (Fig. |l]b-e). A 
coarse graining over the scale of lowest order configurations of this kind is necessary to see 



that the densities of dead and alive perimeter sites are proportional. Reference [19| finds a 
value (0.365 ± 0.01) for the ratio of the number of dead perimeter sites to the total number 
of perimeter sites for 2D square-lattice DLA. 



APPENDIX B: CONSEQUENCES OF NORMALIZATION FOR THE 

MULTIFRACTAL ANALYSIS 

Here, we will briefly explore the consequences of normalizing the growth probabilities Pi 
within an annulus x prior to performing the multifractal analysis suggested in Sec. |V[ We 
determine an x and M dependent normalization factor M{x,M) to multiply each p iy such 
that the sum over the growth probabilities for fixed x equals 1. The normalization procedure 
alters the value of a associated with each pi to 

a { = log pt / log R- log Af(x,M)/ log R. (Bl) 

Now, the distribution of the a is N(a, x, M), which is related to N(a, x, M) by 

N{a,x,M) = N(a - log J\f{x,M)/ log R,x,M). (B2) 

In analogy to Eq. ( pl|) we denote the q th "moment" of N(a,x, M) by Z(q, x,M), i.e., 
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Z(q,x,M) = J daN(a,x,M)M- qu& (B3) 

= wb^! daf,(a ' x ' M)M ' va m 
= wb^w z[q ' x ' M) - (B5) 

From the normalization we know that the first moment of N(a, x, M) is equal to 1. It follows 
that Af(x, M) = Z(l, x, M). 

We continue along the lines in Sec. |V| and define the equivalent f(q,x) to r(q, x), 

~, , v 9 log Z(q,x,M) 

r{q,x) = - lim — B6 

R^oo a log R 

d log Af(x,M) 

= r{q, x)+q hm — . B7 

r-*oo a log R 

If the normalization constant J\f(x, M) displays power-law scaling with the cluster size R, 
then f(q, x) and r(q, x) differ only by a linear function. How does this difference affect 
the Legendre transform of f(q, x) when we compare it to the transform of r(q, x) which is 
defined in Eq. (|55D? First, we calculate the slope oci of r(q, x), 

df(q,x) d\ogAf(x,M) 

a L = — 7^ = a L + hm — . B8 

oq R-^oo d log R 

Since Af(x, M) is always less than one, we see that is merely shifted by a constant to 
smaller values. Second, Legendre transforming f(q, x) yields 

= qa L - f(q, x) (B9) 
d log AT(x,M) 

= qa + q hm — BIO 

H^oo d log R 

dlogX(x,M) 

-T{q, x)-q hm — 

i?->oo Q log R 

= qa -r(q,x) = f L (a L ). (Bll) 

Thus, we retain the functional form of the Legendre transform /l(q;l) of r(q,x) and the 
only difference to /(q;l) is that the latter is shifted towards smaller values of a. 

For example, if we use the formalism presented in this appendix to calculate the /l(q!l) 
for the strong localization case (cf. Sec. |VB|) , we find that 



Jl(&l) = fu and « = h, ( B12 ) 

where denotes the constant value of /^(a(x)) for a given x. In order to interpret the 
result ( pi2|) , consider the normalization condition J2Pi{ x ) — 1 within an annulus. In the 
strong localization case, the annulus x is characterized by only one a. Thus, the the product 
j^-uaj^uf(a) mus t b e constant, leading to a — fz. 
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TABLES 



TABLE I. Summary of multifractal and multiscaling features arising in cases (A), (B) and (C). 







Pa(x,M) 


P(x,M) 


N(a : 


x,M) 






multiscaling 


multiscaling 


multiscaling 


multifractal 


(AH 


= A(a) 


no 


no 


no 


yes 




= A{a)/M» 


yes 


yes 


yes 


no 


(CH 


= A(a)/VlnM 


yes 


yes 


yes 


yes 
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FIGURES 

FIG. 1. (a) The lowest order configuration which contains a dead site (x). (b-e) Examples 
of higher order configurations containing dead sites. 

FIG. 2. (a) Dead and (b) alive growth sites of a DLA cluster with M = 20 000, indicating 
the similar spatial distribution of both types of perimeter sites. Dead sites are perimeter sites 
with pi = 0. Ref. |Rj finds that the number of dead perimeter sites in DLA is proportional to 
the number of all perimeter sites. In our case, dead sites constitute a fraction of « 42% of the 
perimeter sites. 

FIG. 3. Comparison of the two distributions of dead sites, p^\x,M) (broken line), and of 
alive sites, p^(x,M) (solid line). To demonstrate that both are distributed in a similar fashion, 
we have scaled p^(x,M) by the ratio of the number of alive to dead perimeter sites (= 1.36). 

FIG. 4. N(a,x, M) averaged over 18 clusters of mass M = 20 000. Different symbols denote 
different values of a: 1.1(0); F5(D), 1.9(A), 2.3(y) and 2.8(»). To test whether the x dependence 
of N(a, x, M) can be represented by a Gaussian, we plot ±(x — (x)) 2 on the abscissa, where the + 
sign applies if x > (x) and the — sign otherwise. The ordinate scale is logarithmic. Thus, Gaussian 
behavior manifests itself in two straight lines emanating from x — (x) = with slopes of opposite 
sign but equal magnitude. The two solid lines in the plot illustrate this behavior and are intended 
as guides to the eye for the case a = 1.1. 

FIG. 5. Location of a sites from 18 off-lattice DLA clusters of M = 20 000. In (a) we have 
1.5 < a < 1.9, in (b) 2.8 < a < 3.0 and in (c) a > 6. 

FIG. 6. (a) The width £(a, M) and (b) the mean position (x(a,M)) vs M for several values 
of a. Different symbols denote different a values, 0.68 (O), 1-36 (□), 2.04 (A), 2.72 (y), 3.40 (•), 
4.08 (■), 4.76 (filled upward triangle) and 5.95 (filled downward triangle). The data are averaged 
over 18 off-lattice DLA clusters. In (c) (x(a, M)) is plotted as a function of a. Here, different line 
styles correspond to different cluster masses (see legend). 
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FIG. 7. Dependence on q of r(q,x) for different values of x. Here, r(q,x) was determined 
by fitting a straight line through the 3 data points corresponding to M = 5 000, 10 000, 20 000 
in a plot of log Z(q, x, M) vs logM averaged over 18 DLA clusters. Different line styles denote 
the different x values (legend). For comparison, we plot the r(q) resulting from an analysis of the 



growth probabilities of the entire cluster (O)- Since 2D DLA displays a phase transition HQ, 41 



r(q, x) for negative q can only be considered an "effective" exponent which will display larger and 
larger slopes as M — > oo. 
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